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Abstract:  The  hyperbolic  and  plasmonic  properties  of  silicon  nanowire/ Ag 
arrays  have  been  investigated.  The  aligned  nanowire  arrays  were  formed 
and  coated  by  atomic  layer  deposition  of  Ag,  which  itself  is  a  metamaterial 
due  to  its  unique  mosaic  film  structure.  The  theoretical  and  numerical 
studies  suggest  that  the  fabricated  arrays  have  hyperbolic  dispersion  in  the 
visible  and  IR  ranges  of  the  spectrum.  The  theoretical  predictions  have  been 
indirectly  confirmed  by  polarized  reflection  spectra,  showing  reduction  of 
the  reflection  in  p  polarization  in  comparison  to  that  in  s  polarization. 
Studies  of  dye  emission  on  top  of  Si/Ag  nanowire  arrays  show  strong 
emission  quenching  and  shortening  of  dye  emission  kinetics.  This  behavior 
is  also  consistent  with  the  predictions  for  hyperbolic  media.  The  measured 
SERS  signals  were  enhanced  by  almost  an  order  of  magnitude  for  closely 
packed  and  aligned  nanowires,  compared  to  random  nanowire  composites. 

These  results  agree  with  electric  field  simulations  of  these  array  structures. 
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1.  Introduction 

Metamaterials  are  engineered  composite  materials  containing  subwavelength  inclusions  with 
rationally  designed  shapes,  sizes,  mutual  arrangements  and  orientations,  which  often  exhibit 
unparalleled  properties  not  available  in  nature  or  traditional  synthetic  materials  [1,2].  The 
field  of  metamaterials  has  brought  about  novel  designs,  since  researchers  have  been  able  to 
overcome  the  constraints  associated  with  conventional  materials  and  explore  entirely  new 
concepts  with  a  much  larger  range  of  material  parameters  available  through  structured  media. 
Materials  with  negative  indices  of  refraction  [3-7],  “perfect”  lenses  [8,9]  and  hyperlenses 
[10-14]  that  focus  and  image  with  resolution  beyond  the  diffraction  limit,  invisibility  cloaks 
[15,16],  and  nanoscopic  lasers  [17-21]  are  just  some  of  the  novel  concepts  that  have  been  put 
forward. 

A  new  class  of  uniaxial  (meta)materials  has  been  reported  [22-25],  which  have  highly 
unusual  iso-frequency  surfaces  defined  by  allowed  wavevectors  at  constant  frequency.  When 
all  principal  components  of  the  dielectric  permittivity  tensor  are  positive,  the  iso-frequency 
surface  is  “closed”  and  forms  a  spheroid  or  ellipsoid  in  the  wavevector  space  (k- space).  In 
such  media,  the  magnitude  of  the  wavenumber  is  limited,  constraining  the  resolution  of 
optical  imaging,  and  the  density  of  photonic  states.  In  contrast,  a  highly  anisotropic 
metamaterial,  whose  dielectric  permittivity  in  orthogonal  directions  has  different  signs, 
exhibits  a  hyperbolic  iso-frequency  surface,  and  thus  is  referred  to  as  hyperbolic.  Such 
materials  have  an  almost  unlimited  magnitude  of  a  wave  vector  k  and  near  infinite  local 
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density  of  photonic  states  [26].  Recent  theoretical  predictions  [26]  and  experimental  results  in 
hyperbolic  metamaterial  structures  [27-33]  have  shown  that,  in  line  with  the  seminal  work  by 
Purcell  [34],  the  high  density  of  photonic  states  can  control  spontaneous  emission  of 
luminescent  centers  and  reflectance  of  roughened  metamaterials’  surfaces  [35,36]. 

The  most  common  approach  to  fabricate  a  metamaterial  is  by  patterning  periodic  arrays  of 
individual  elements,  which  can  be  metal,  dielectric,  plastics,  etc.,  depending  on  application. 
The  properties  of  metamaterials  can  then  be  tailored  by  controlling  the  individual  element 
shape,  size  and  geometry,  as  well  as  the  array  periodicity  and/or  arrangement  [1,2].  However, 
most  of  traditional  cleanroom  fabrication  techniques  are  prohibitively  expensive  and  hardly 
scalable  [37].  Much  more  appealing  are  bottom-up  techniques  capable  of  producing 
centimeter  size  samples,  which  include  electroplating  growth  of  metallic  nanowires  in  porous 
alumina  membranes  [38-42],  eutectic  crystal  growth  of  split  ring  resonator  structures  [43], 
and  deposition  of  lamellar  semiconductor  or  metal-dielectric  film  structures  [28-32,44].  In 
addition,  a  metamaterial  can  also  be  formed  by  Atomic  Layer  Deposition  (ALD)  growth- 
induced  patterning,  as  has  recently  been  reported  for  Ag  [45]. 

The  density  of  photonic  states  in  hyperbolic  metamaterials  (which  can  control  a  variety  of 
quantum  and  classical  phenomena)  is  significantly  larger  inside  the  structure  than  on  its 
surface  [29].  Therefore,  many  applications  require  impregnation  of  hyperbolic  metamaterials 
with  nonlinear  and  luminescent  chromophores  and  other  active  components.  One  realization 
of  hyperbolic  metamaterials,  allowing  their  easy  functionalization,  is  arrays  of  parallel  metal- 
coated  nanowires  that  are  separated  by  air  gaps,  which  is  studied  in  this  work  theoretically 
and  experimentally. 

In  this  work,  we  fabricate  and  examine  the  optical  properties  of  closely  spaced  long  Si/Ag 
nano  wire  composite  arrays.  We  predict,  theoretically  and  numerically,  that  these  structures 
have  hyperbolic  dispersion  in  the  visible  and  infrared  ranges  of  the  spectrum.  The  results  of 
reflection  and  dye  emission  experiments  are  consistent  with  the  model  predictions.  We  also 
report  an  enhancement  of  surface  enhanced  Raman  scattering  (SERS),  which  is,  most  likely, 
due  to  plasmonic  coupling  of  closely  spaced  nanowires  and  enhanced  density  of  photonic 
states  in  hyperbolic  metamaterials. 

2.  Formation  of  closely  spaced  Si  nanowire  arrays 

There  are  a  number  of  approaches  in  forming  silicon  nanowires  (NWs),  such  as  chemical 
vapor  deposition  [46]  and  thermal  evaporation  [47],  but  in  general,  these  processes  require 
very  specific  growth  chambers  which  are  expensive  and  which  can  require  rather  dangerous 
gases,  such  as  silane.  Furthermore,  most  of  these  systems  produce  randomly  grown  Si 
nanowires.  In  order  to  form  aligned  and  ordered  Si  nanowire  arrays  with  small  diameters,  the 
vapor  liquid  solid  (VLS)  growth  mechanism  is  used  [48],  and  small  diameter  nanowires  can 
be  formed  by  control  of  the  metal  catalyst  size,  which  may  require  further  lithographic  steps. 
Instead,  we  have  chosen  a  metal  assisted  chemical  etching  process  [49],  which  allows  a 
simple  and  inexpensive  way  of  forming  silicon  nanowire  arrays,  with  very  few  equipment 
requirements. 

Preparation  of  silicon  nanowires  through  chemical  etching  has  been  found  to  work  well 
when  silver  is  used  to  enhance  the  etching  [50].  The  etching  process  can  be  described  as 
follows:  Ag  nanoparticles  are  deposited  onto  the  silicon  substrate.  Some  of  the  Ag  is  oxidized 
and  dissolved  in  solution,  forming  Ag  +  ions,  which  react  with  Si  at  the  Si/Ag  nanoparticle 
interface,  so  that  enhanced  silicon  etching  occurs,  and  it  is  localized  to  the  regions  covered  by 
the  Ag  nanoparticles.  As  this  enhanced  Si  etching  continues,  channels  are  formed  into  the 
silicon  substrate,  leaving  silicon  walls  in  regions  with  no  Ag  nanoparticles,  thus  creating  a  Si 
NW  array. 

A  number  of  Silicon  substrates  were  investigated,  including  Si(100),  Si(lll),  lightly 
doped  p  and  n-type  wafers,  heavily  doped  p  and  n-type  wafers,  as  well  as  p-doped  50  micron 
thick  Si  wafers.  The  silicon  wafers  were  cut  into  lxl  cm2  pieces,  washed  with  acetone, 
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ethanol,  and  DI  water,  followed  by  a  piranha  etch  (sulfuric  acid  and  hydrogen  peroxide  in  a 
3:1  v/v  ratio).  The  samples  were  then  dipped  in  5%  HF/water  solution  for  3  minutes,  placed  in 
a  4.8  M  HF  and  0.005M  AgN03  solution  to  form  the  Ag  catalyst  [49].  The  resultant  samples 
were  rinsed  in  deionized  (DI)  water  and  the  remaining  Ag  catalyst  was  removed  using  a  nitric 
acid/DI  solution.  The  nano  wires  were  characterized  via  scanning  electron  microscopy  (SEM). 

2. 1  Effect  of  Si  substrate  doping  on  array  formation 

A  number  of  silicon  substrates  were  investigated,  using  the  optimal  electroless  Ag  catalyst 
deposition  and  optimal  HF  etching  conditions.  Figure  1  shows  the  SEM  images  of  some  of 
the  resulting  nanowires  etched  under  the  same  conditions  as  a  function  of  substrate  doping. 
Figure  1(a)  shows  an  SEM  micrograph  of  a  lightly  boron  doped  p-type  silicon  wafer  with  a 
resistivity  of  1-30  ohm-cm,  Fig.  1(b)  shows  an  image  for  a  lightly  phosphorus  doped  n-type 
silicon  with  a  resistivity  of  1-20  ohm-cm,  Fig.  1(c)  shows  a  heavily  boron  doped  p-type 
silicon  wafers  with  resistivity  of  0.01-0.02  ohm-cm,  Fig.  1(d)  shows  a  heavily  arsenic  doped 
n-type  silicon  with  a  resistivity  of  0.0015-0.007  ohm-cm  and  Fig.  1(e)  shows  Si  arrays  that 
can  be  formed  on  both  sides  of  a  Si  wafer,  assuming  the  etching  conditions  are  properly 
optimized. 


Fig.  1.  Si  NWs  produced  using  a)  p-type  silicon,  1-30  ohm-cm,  b)  n-type  silicon,  1-20  ohm- 
cm,  c)  p-type  silicon,  0.01-0.02  ohm-cm,  d)  n-type  silicon,  0.0015-0.007  ohm-cm  and  e)  p-type 
Si  (100)  50  micron  thick  wafer,  etched  from  both  sides. 


As  is  evident  from  Fig.  1,  the  metal  assisted  chemical  etching  process  worked  well  for 
most  types  of  silicon  substrates  except  those  that  were  heavily  doped.  This  is  evident  in  Fig. 
1(d)  for  a  very  heavily  doped  sample,  which  resulted  in  short,  disordered  nanowires  that  were 
not  properly  separated  from  each  other. 

3.  Formation  of  NW/Ag  composite  arrays 

In  order  to  investigate  both  hyperbolic  and  plasmonic  properties  of  the  NW  arrays,  we  coated 
the  Si  NW  arrays  with  21  nm  of  Ag,  using  atomic  layer  deposition  (ALD).  The  ALD  process, 
due  to  its  conformality  and  high  aspect  ratio  coverage,  was  deemed  ideal  for  achieving  full 
coverage  of  the  closely  spaced  NW  arrays. 

The  growth  of  the  Ag  films  was  performed  by  remote  plasma  enhanced  atomic  layer 
deposition  (PEALD),  using  Ag(fod)(PEt3)  (fod  =  2,2-dimethyl-6,6,7,7,8,8,8- 
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heptafluorooctade-3,5-dionato)  and  hydrogen  plasma  as  the  precursors  [45,  51].  The  Ag  films 
were  deposited  at  125°C  onto  the  Si  nanowire  arrays,  as  well  as  on  random  Si  nano  wires. 


Fig.  2.  a)  SEM  image  of  21  nm  thick  Ag  thin  film  deposited  on  Si  using  plasma  enhanced 
atomic  layer  deposition  (PEALD);  b)  SEM  image  of  Si  NW  arrays  coated  with  ALD  Ag  near 
the  substrate. 

The  thicknesses  of  the  deposited  Ag  films  were  7  nm,  17  nm,  21  nm  and  38  nm,  as 
measured  by  ellipsometry  and  XTEM. 

The  PEALD  film  consisted  of  poly  crystalline  fee  flat  Ag  islands  (Fig.  2(a)),  separated  by 
very  small  air  gaps.  No  detectable  contaminants  (from  EDS)  were  noted  and  the  thicker  films 
had  a  resistivity  in  line  with  metallic  Ag.  In  terms  of  coating  the  Si  NW  arrays,  although  the 
Si  nano  wires  in  the  arrays  were  30  microns  long  and  relatively  closely  spaced  (separated  by 
only  150  -  200  nm),  we  see  conformal  coverage  of  the  ALD  Ag  along  the  whole  length,  down 
to  the  substrate,  as  shown  in  Fig.  2(b). 

In  terms  of  the  optical  response,  however,  the  ALD  Ag  does  not  resemble  a  standard  Ag 
thin  film,  in  that  the  ALD  Ag  exhibits  plasmonic  behavior,  as  reported  earlier  [45].  This 
plasmonic  behavior  is  also  evident  from  Fig.  3(a),  which  shows  an  experimentally  determined 
reflection  spectrum  for  the  PEALD  Ag,  along  with  an  FDTD  simulation  of  this  Ag  mosaic 
structure.  The  resonance  near  600  nm  is  due  to  plasmonics,  and  it  has  been  shown  to  blue 
shifts  with  decreasing  film  thickness  [52].  This  resonance  is  only  present  in  ALD  Ag  and  not 
in  Ag  films  deposited  by  other  means  [52],  which  agrees  well  with  the  SERS  enhancement 
only  seen  for  ALD  Ag  [45].  Thus  PEALD  Ag  is  in  fact  a  metamaterial  due  to  its  unique 
mosaic  grain  structure,  modeled  as  a  coaxial  cylinder  of  Ag  surrounded  by  an  air  gap  [45]. 
Figure  3(b)  shows  a  single  coaxial  cylinder  used  to  model  the  Ag  ALD  film  mosaic 
morphology  in  the  FDTD  simulation.  The  simulation  was  obtained  by  forming  an  array  of 
these  cylinders 


wavelength  (nm) 


Coaxial  cylinder 


Fig.  3.  a)  Comparison  of  the  measured  reflection  data  of  ALD  Ag  to  a  finite  difference  time 
domain  simulation  of  a  200nm  diameter  air  ring  in  a  40nm  thick  Ag  film.  Note  that  the  data 
and  simulation  both  have  the  same  line  shape  and  that  the  positions  of  the  resonance  agree 
well,  b)  Electric  field  simulation  of  a  single  element  of  the  PEALD  Ag  microstructure  (coaxial 
cylinder)  used  to  obtain  the  simulation  in  a).  Also  shown  are  the  resultant  high  electric  fields  in 
and  above  the  air  gap. 
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in  order  to  obtain  the  reflection  data  shown  in  Fig.  3(a)  and  the  electric  field  enhancements 
shown  in  red  in  Fib  3(b).  Note  that  these  high  electric  fields  form  not  only  in  the  air  gaps  but 
also  above  the  gaps,  which  leads  to  strong  plasmonic  behavior  as  measured  by  SERS  [45]. 

4.  Theoretical  prediction  of  Hyperbolic  Dispersion 

To  understand  the  topological  properties  of  the  waves  propagating  in  our  metamaterials,  we 
used  2D  fmite-element-based  solutions  of  Maxwell  equations  [53]  to  deduce  the  dependence 
of  the  z-component  of  wavevector  on  frequency  and  x-component  of  the  wavevector.  In  these 
calculations  we  used  a  more  general  model  of  the  structure  than  the  one  realized  in 
experiments.  Thus,  the  system  was  represented  as  a  periodic  array  of  infinitely  long 
composite  wires  with  shelled  internal  structure  (one  unit  cell  of  such  structure  is  shown  in 
Fig.  4(a)).  These  calculations  have  revealed  that,  similar  to  previously  studied  nanowire-based 
metamaterials  [38-42],  Si-based  nanowire  arrays  support  two  waves  that  differ  by  their 
polarization.  The  dispersion  of  these  two  waves  was  numerically  calculated  for  the  following 
set  of  parameters:  free-space  wavelength  900  nm,  distance  between  centers  of  Si  nanowires 
250  nm,  the  diameter  of  nanowires:  150  nm,  with  thickness  of  ALD  silver  coatings  ranging 
from  7  nm  to  38  nm.  In  these  calculations,  we  assumed  that  optical  properties  of  components 
of  metamaterials  are  significantly  close  to  those  of  bulk  media  [54]. 

The  calculations  show  that  for  small  angles  of  incidence,  the  behavior  of  the  metamaterial 
is  consistent  with  what  is  expected  of  a  uniaxial  homogeneous  medium.  For  smaller  values  of 
silver  thickness,  the  material  exhibits  elliptical  behavior,  while  for  larger  thickness  of  the 
silver  shells,  the  behavior  is  hyperbolic-like.  Least-square  fitting  of  these  numerically 
calculated  dispersion  curves  to  ellipses  or  hyperbolae  yields  the  effective  medium  parameters 
for  the  structure  [40].  These  effective  medium  parameters  are  summarized  in  Fig.  4(b). 

In  the  limit  when  the  nanowire  radius,  shell  thickness,  and  distance  between  the  wires  are 
much  smaller  than  the  wavelength,  the  dispersion  of  nanowire  structures  can  be  described  by 


Fig.  4.  (a)  schematic  geometry  of  the  unit  cell  for  a  multi-shell  nano  wire  composite,  used  in 
derivation  of  effective  medium  response;  experimental  configuration  corresponds  to  a  two 
shell  system  or  N  =  2;  (b)  effective  medium  parameters  extracted  from  finite-element  solutions 
of  Maxwell  equations  (symbols)  and  from  effective  medium  theory  (lines);  (c)  spectral 
dependence  of  effective  medium  parameters  for  different  thickness  of  Ag  shell;  (d)  dependence 
of  effective  plasma  wavelength  of  composite  in  (c)  as  a  function  of  Ag  shell  thickness. 

effective  medium  theory  that  can  be  calculated  by  straightforward  generalization  of  Maxwell- 
Gamett  approach  to  describe  the  composite  containing  N  shell-type  components.  Explicitly, 
the  component  of  the  permittivity  tensor  along  the  wire,  8ZZ,  is  calculated  as  the  weighed 
average  of  the  permittivities  of  the  constituent  components  as 
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with  pi  being  the  (surface)  concentration  of  the  i-th  component  of  the  composite. 

To  calculate  the  component  of  permittivity  tensor  in  the  direction  perpendicular  to  the 
wires,  it  is  necessary  to  find  field  distribution  along  the  nanowire  system.  Straightforward 
calculation  yields: 
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The  agreement  of  the  numerically-extracted  effective  medium  parameters  with  prediction 
of  generalized  Maxwell- Garnett  theory  is  illustrated  in  Fig.  4(b).  It  is  seen  that  the  iterative 
technique  adequately  describes  the  behavior  of  e±  in  the  regime  when  the  concentration  of 
inclusions  is  relatively  small  (the  relative  concentration  of  the  host  component  of 
metamaterial  pN  >70%).  In  contrast,  the  parallel  component  of  the  permittivity  is  only 
qualitatively  described  by  the  effective  medium  theory.  This  apparent  disagreement  has  its 
origin  in  strong  spatial  dispersion  that  is  expected  in  the  vicinity  of  epsilon-near-zero 
behavior  for  composites  with  relatively  large  unit  cells.  Similar  deviations  have  been 
previously  observed  in  gold  nano  wire  metamaterials  [55].  The  implications  of  this  behavior 
for  optical  response  of  the  nanowire  composites  and  the  analytical  description  of  the 
properties  of  nano  wire  systems  in  ENZ  regime  will  be  detailed  in  our  future  works. 

To  explore  the  parameter  space  enabled  by  the  Si/Ag  nano  wire  arrays,  effective  medium 
theory  was  used  to  analyze  the  dependence  of  effective  medium  parameters  as  a  function  of 
wavelength  for  several  fixed  values  of  the  Ag  shell  thickness.  The  results  of  these  calculations 
are  summarized  in  Fig.  4(c).  It  is  seen  that,  similar  to  conventional  plasmonic  nanowire 
arrays,  Si/Ag  composites  offer  broadband  hyperbolic  response,  with  plasma-like  behavior  of 
By.  As  expected,  the  thickness  of  plasmonic  shell  plays  a  role  of  a  convenient  control 
parameter,  that  can  be  used  to  tune  effective  plasma  frequency  (the  frequency  at  which  By 
crosses  zero).  Such  tuning  is  shown  in  Fig.  4(d). 


5.  Experimental  Studies  of  Reflection  and  Emission  in  Si/Ag  arrays 

5. 1  Specular  and  Diffused  Reflection  of  Si/Ag  arrays 

The  four  nanowire  samples  with  the  parameters  described  in  Section  4  and  thickness  of  ALD 
silver  equal  to  7  nm,  17  nm,  21  nm  and  38  nm  have  been  studied  experimentally.  According 
to  the  theoretical  modeling  presented  in  Section  4,  all  these  samples  have  hyperbolic 
dispersion  in  the  red  and  infrared  ranges  of  the  spectrum.  All  of  them  looked  dark  gray,  likely 
due  to  a  combination  of  roughness,  which  is  inherent  to  nanomaterials  fabricated  by  etching, 
and  (as  explained  below)  preferential  scattering  of  light  inside  roughened  hyperbolic 
metamaterials  [35,  36].  Correspondingly,  the  specular  reflection  of  the  samples  was  very  low, 
under  0.5%.  We  have  attempted  to  obtain  the  values  of  dielectric  permittivities  extracted  from 


#187087  -  $15.00  USD  Received  15  Mar  2013;  revised  15  May  2013;  accepted  16  May  2013;  published  17  Jun  2013 
(C)  2013  OSA  17  June  2013  |  Vol.  21,  No.  12  |  DOI:10.1364/OE.21.014962  |  OPTICS  EXPRESS  14969 


the  angular  reflectance  data  (following  [42]).  However,  due  to  the  very  weak  specular 
reflection  and  strong  diffused  scattering,  it  is  not  clear  if  meaningful  results  can  be  obtained 
by  this  method. 


Fig.  5.  Reflection  spectra  of  Si/ALD  silver  sample  with  21  nm  ALD  silver  coating,  measured 
in  p  and  s  polarizations  in  an  integrating  sphere  setup  at  -10°  incidence  angle. 

At  the  same  time,  the  reflection  measured  in  an  integrated  sphere  setup  (which  captured 
both  specularly  reflected  light  and  diffusely  reflected  light)  was  significantly  larger;  see  for 
example  Fig.  5,  depicting  the  reflection  spectrum  of  the  sample  with  21  nm  ALD  silver 
coating.  One  can  see  that  the  reflectivity  in  p  polarization  is  significantly  smaller  than  that  in  s 
polarization,  in  particularly  in  the  near-infrared  part  of  the  spectrum.  This  is  expected  from 
rough  surfaces  of  metamaterials  with  hyperbolic  dispersion  [35,36],  in  which  enhanced 
scattering  of  p  polarized  light  inside  the  hyperbolic  medium  (where  the  density  of  photonic 
states  is  high)  causes  substantial  reduction  of  reflection.  No  substantial  reduction  of 
reflectance  is  predicted  for  s  polarized  light  [36].  Correspondingly,  the  experimental  result 
depicted  in  Fig.  5  is  consistent  with  the  theoretical  prediction  of  hyperbolic  dispersion  in  the 
sample. 

5.2  Kinetics  and  intensity  of  dye  emission  on  top  of  Si/Ag  nanowire  array 

In  this  particular  experiment,  we  studied  the  array  of  Si/Ag  nano  wires  coated  with  21  nm 
ALD  silver  layer.  A  film  of  polymethyl  methacrylate  (PMMA)  doped  with  IR140  laser  dye 
(in  concentration  0.013  M)  was  deposited  on  top  of  the  nanowire  array  sample  and  onto  the 
control  substrates,  which  were  glass  and  silver  film  deposited  on  glass.  The  film  thickness  in 
control  samples  was  -80  nm.  The  Si/Ag  nano  wire  sample  had  approximately  the  same 
concentration  of  dye  molecules  per  unit  area  of  the  sample.  However,  dye  penetrated  in  the 
voids  between  nanowires,  and  the  effective  thickness  could  not  be  determined. 

In  order  to  measure  the  emission  kinetics,  the  IR140:PMMA  films  were  pumped  using  a 
100  fs  Ti:  sapphire  laser  at  X  =  792  nm.  The  emission  was  detected  using  a  near- infrared 
Hamamatsu  streak  camera  equipped  with  appropriate  long-pass  filters  (2  >850  nm).  All 
emission  kinetics,  averaged  over  multiple  measurements,  along  with  the  kinetics  of  the 
pumping  pulse  scattered  by  a  pure  glass  substrate  (which  shows  the  time  resolution  of  the 
apparatus  determined  by  jitter  of  the  laser  and  wide  open  slit  of  the  streak  camera)  are 
depicted  in  Fig.  6.  One  can  see  that  the  emission  of  dye  on  top  of  the  metamaterial  sample  is 
strongly  quenched  in  first  -100  ps,  and  its  decay  kinetics  is  significantly  shortened  (more  than 
twofold  in  comparison  to  the  control  samples)  at  longer  times. 
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Fig.  6.  (a)  Spontaneous  emission  kinetics  of  the  IR140:PMMA  film  deposited  on  the  top  of 
glass  (1),  200  nm  thick  silver  film  (2),  and  the  Si/ALD  Ag  metamaterial  sample  (3).  Trace  4 
shows  the  time  resolution  of  the  apparatus. 

We  have  also  studied  the  emission  spectra  of  IR140:PMMA  films  deposited  on  top  of  the 
metamaterial  sample  and  the  glass  substrate.  In  this  particular  experiment,  the  samples  were 
pumped  at  784  nm.  The  emission  signal  was  detected  with  a  photomultiplier  tube  (PMT) 
connected  to  the  exit  slit  of  a  monochromator.  The  emission  spectra  are  shown  in  Fig.  7.  It 
can  be  seen  that  the  emission  of  IR140:PMMA  deposited  on  top  of  the  metamaterial,  is 
almost  one  order  of  magnitude  smaller  compared  to  the  same  dye-doped  film  deposited  on 
glass.  (The  difference  can  be  even  larger  considering  the  large  error  bar.) 


Fig.  7.  Spontaneous  emission  spectra  of  the  IR140:PMMA  film  deposited  on  top  of  the  Si/Ag 
metamaterial  sample  with  21  nm  ALD  silver  coating  (1)  and  glass  (2),  pumped  at  X  =  784  nm 
into  the  absorption  band  of  IR140. 

The  experimentally  observed  reduction  of  the  dye  emission  intensity  and  shortening  of  the 
emission  kinetics  are  consistent  with  the  theoretical  predictions  for  hyperbolic  metamaterials 
[26].  In  fact,  high  density  of  photonic  states  on  the  surface  (and  inside)  of  a  metamaterial  with 
hyperbolic  dispersion  causes  enhancement  of  spontaneous  emission  directed  inside  the 
material.  This  determines  the  emission  kinetics  shortening  and  reduction  of  the  emission 
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intensity.  Similar  emission  quenching  and  shortening  of  the  emission  decay  time  have  been 
observed  in  a  dye-doped  film  deposited  on  top  of  array  of  silver  nanowires  grown  in  alumina 
membrane  [27]. 

5.3  Angular  distribution  of  dye  emission  in  Si/Ag  nanowire  array 

Hyperbolic  metamaterials  have  been  predicted  to  produce  cone-shaped  or  beam-like  patterns 
of  spontaneous  emission  inside  the  medium  [26].  The  angular  distribution  of  emission  outside 
the  sample  is  much  less  researched.  In  the  experimental  study  of  this  phenomenon,  we 
deposited  PMMA  film  doped  with  HITC  laser  dye  (in  concentration  0.04  M)  onto  array  of 
Si/Ag  nanowires  with  38  nm  silver  coating  as  well  as  onto  ALD  silver  film  and  thermal  vapor 
deposited  silver  film.  The  thickness  of  the  dye-doped  PMMA  films  on  top  of  silver  films  was 
~80  nm.  In  the  Si/Ag  nanowire  sample,  dye  (which  had  approximately  the  same  concentration 
of  dye  molecules  per  unit  area)  penetrated  between  the  wires. 

Dye  molecules  were  excited  with  -5  nm  pulses  of  the  optical  parametric  oscillator  (OPO) 
at  X  =  750  nm  and  45°  incidence  angle.  The  emission  intensity  was  measured,  in  the  spectral 
band  centered  at  815  nm,  at  a  variety  of  detection  angles  ranging  from  0°  (normal  to  the 
sample  surface)  to  -90°  (parallel  to  the  surface).  A  long-pass  filter  in  front  of  the  detector 
blocked  laser  pumping  and  transmitted  dye  emission. 

Angular  distributions  of  emission  in  three  samples  studied  are  depicted  in  Fig.  8.  One  can 
see  that  the  emission  patterns  from  the  Si/Ag  nanowire  array  and  ALD  silver  film  are  very 
close  to  Lambertian  (~cos(0)).  This  behavior  could  probably  be  expected,  since  surfaces  of 
both  samples  are  rough  (the  roughness  of  Si/Ag  NW  array  much  larger  than  that  of  ALD 
silver.)  At  the  same  time,  the  emission  of  dye  on  top  of  a  much  smoother  (>5nm)  thermal 
vapor  deposited  silver  is  slightly  broader.  The  comparison  of  the  experiment  with  the  theory 
will  be  published  elsewhere. 


Angle  (degree) 


Fig.  8.  Angular  distributions  of  emission  of  HITC  dye  molecules  on  top  of  Si/Ag  NW  array 
(diamonds),  ALD  silver  film  (triangles)  and  silver  film  deposited  via  thermal  vapor  deposition 
(squares).  Solid  line:  cos(0). 

6.  Plasmonic  Properties  of  Nanowire  Arrays 

Raman  scattering  is  a  technique,  which  can  measure  interatomic  vibrations  that  are  unique  to 
each  material,  thus  providing  a  chemical  “fingerprint”  for  any  material.  Unfortunately,  the 
scattering  efficiencies  are  very  low,  roughly  10”8  of  the  incident  laser  light.  However,  the 
signal  can  be  greatly  enhanced  by  the  use  of  metal  nanoparticles,  which  enhance  the  local 
electric  fields  [56].  This  process  is  referred  to  as  Surface  Enhanced  Raman  Scattering 
(SERS).  One  type  of  SERS  nanostructures  resulting  in  strong  SERS  signal  are  dielectric 
core/metal  sheath  nanowires  [57,58],  which  have  been  shown  to  create  electric  field  hotspots 
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at  their  crossings  due  to  plasmonic  coupling.  It  has  also  been  shown  that  two  closely  spaced 
parallel  dielectric  core/metal  sheath  nano  wires  result  in  a  greater  enhancement  than  those 
crossed  at  90°  [58]. 

A  self-assembled  monolayer  of  benzenethiol  (BZT)  was  deposited  on  the  NW  arrays  by 
soaking  them  in  10~4  M  BZT/toluene  solution  overnight  and  rinsing  with  methanol,  which 
removed  excess  BZT  not  bonded  to  the  Ag  surface.  The  SERS  spectra  were  obtained  using  a 
DeltaNu  Reporter  785  micro-Raman  system,  operating  at  785  nm  and  the  spectra  were 
collected  with  a  laser  power  of  less  than  3  mW.  In  addition,  COMSOL  electric  field 
calculations  were  performed  in  order  to  interpret  the  experimental  SERS  results.  The  SERS 
results  for  the  Si/Ag  NW  arrays  are  shown  in  Fig.  9(a).  The  major  Raman  peaks  can  be 
assigned  to  symmetric  ring  breathing,  in  plane  C-H  bending  and  in  plane  C-C  stretching  of 
the  phenyl  ring  from  the  active  molecules,  benzenethiol  (BT),  which  are  in  good  agreement 
with  those  reported  previously  [59,  60].  As  can  be  seen,  very  strong  lines  associated  with  the 
BZT  molecule  are  seen,  and  the  uniformity  of  the  SERS  signal  intensity  is  quite  uniform 


Fig.  9.  (a)  SERS  spectrum  for  benzenethiol  (BZT)  for  Si/Ag  NW  arrays  compared  to 
horizontal  random  Si/Ag  NWs,  and  (b)  SERS  spectrum  for  three  random  locations  on  the  Si 
NW  array. 

when  taken  from  different  regions  of  the  sample  (Fig.  9(b)).  More  importantly,  the  SERS 
signal  of  the  closely  spaced  Si/Ag  NW  arrays  is  a  factor  of  7  stronger  compared  to  a  similar 
density  of  random  Si/Ag  nanowires  placed  on  a  flat  Si  substrate,  Fig.  9(a).  To  arrive  at  this 
number,  the  SERS  intensity  of  the  arrays  was  normalized  to  the  total  number  of  NWs  in  a  2 
micron  diameter  area  (laser  spot  size).  Similarly,  the  SERS  intensity  of  the  random  horizontal 
nano  wires  was  normalized  to  the  total  number  of  random  NWs  in  the  2  micron  area. 

To  explain  the  enhancement  in  the  SERS  intensity  for  the  closely  spaced  NW  arrays,  we 
refer  to  Kotmann  and  Martin  [61],  who  have  shown  that  2  parallel  NWs  with  a  small 
separation  distances  can  result  in  the  formation  of  hot  spots  corresponding  to  very  intense 
electric  field,  leading  to  a  significantly  enhanced  coupled  plasmon  resonance. 

In  order  to  understand  the  enhancement  in  the  case  of  our  closely  spaced  aligned  high 
density  nanowire  arrays,  we  have  performed  COMSOL  electric  field  calculations,  shown  in 
Fig.  10,  which  show  the  electric  fields  induced  between  the  NW  arrays  upon  irradiation  by 
light.  As  can  be  seen  in  Fig.  10(a),  for  a  4x4  array,  parallel  and  closely  packed  NWs  induced 
high  electric  fields  along  the  whole  lengths  of  the  nanowires,  which  may  account  for  the  very 
large  enhancement  of  the  observed  SERS  signals.  Figure  10(b)  is  a  side  view  of  the  parallel 
NW  arrays  and  the  electric  fields  generated,  clearly  demonstrating  the  strong  plasmonic 
coupling  between  the  NWs  in  the  array.  Due  to  the  difficulty  of  computing  electric  fields  for 
long  NWs,  only  300  nm  long  wires  were  modeled.  However,  these  results  are  also  applicable 
to  longer  NWs.  Furthermore,  COMSOL  simulations  also  show  that 
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Fig.  10.  COMSOL  simulation  of  the  electric  field  enhancement  of  closely  spaced  Ag  nanowire 
arrays  a)  looking  from  top  and  b)  along  the  NWs.  In  the  simulations,  the  NWs  were  100  nm  in 
diameter  and  300  nm  in  length. 

on  average,  the  strong  plasmonic  coupling  is  not  dramatically  influenced  by  small  variations 
in  the  nano  wire  diameters.  Since  the  nanowires  in  the  arrays  are  aligned  but  randomly  spaced 
from  each  other,  if  one  assumes  that  some  fraction  of  the  nanowires  are  significantly  closer 
than  the  average  distance  of  150-200  nm,  strong  plasmon  coupling  and  thus  a  high  SERS 
intensity  would  be  expected,  as  discussed  by  Kotmann  et  al.  [61].  The  fact  that  these 
nano  wires  form  regions  where  they  are  much  more  closely  spaced  is  clearly  seen  in  Fig.  1(b), 
1(c). 

Conclusions 

In  summary,  we  have  synthesized  arrays  of  silver  coated  silicon  nanowires  and  studied  their 
properties  theoretically  and  experimentally.  The  arrays,  which  had  a  thicker  Ag  coating,  have 
been  predicted  to  have  a  hyperbolic  dispersion  in  the  visible  and  infrared  ranges  of  the 
spectrum,  modeled  by  finite-element-based  solutions  of  Maxwell’s  equations.  In  terms  of 
experimental  results,  we  have  noted  a  significant  reduction  of  the  samples’  reflectance  at  p 
polarization  (in  comparison  to  s  polarization),  which  is  expected  from  rough  surfaces  of  a 
hyperbolic  metamaterial.  Furthermore,  the  emission  behavior  of  dye  molecules  deposited  on 
top  of  the  nanowire  arrays  has  been  obtained,  and  the  results  show  strong  dye  emission 
quenching,  by  nearly  an  order  of  magnitude,  and  shortening  of  the  emission  kinetics.  These 
two  observations  are  consistent  with  the  theoretical  predictions  for  hyperbolic  metamaterials. 

The  measured  SERS  signals  for  the  Si/Ag  array  composites  were  enhanced  by  almost  an 
order  of  magnitude  for  closely  packed  and  aligned  nanowires,  compared  to  random  nanowire 
composites  deposited  on  a  flat  surface.  COMSOF  electric  field  simulations  suggest  that  this 
enhancement  is  due  to  plasmon  coupling  in  closely  spaced  nanowires,  which  are  present  in 
the  arrays  investigated. 
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